Static charged perfect fluid spheres in general relativity 
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I. INTRODUCTION 
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Tzarigradsko Shausse 72, Sofia 1784, Bulgaria 

Interior perfect fluid solutions for the Reissner-Nordstrom metric are studied on the basis of a new 
classification scheme. It specifies which two of the characteristics of the fluid are given functions and 
accordingly picks up one of the three main field equations, the other two being universal. General 
formulae are found for charged de Sitter solutions, the case of constant energy component of the 
energy-momentum tensor, the case of known pressure (including charged dust) and the case of 
linear equation of state. Explicit new global solutions, mainly in elementary functions, are given as 
illustrations. Known solutions are briefly reviewed and corrected. 
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The unique exterior metric for a spherically symmetric charged distribution of matter is the Reissner-Nordstrom 
solution. Interior regular charged perfect fluid solutions are far from unique and have been studied by different authors. 
The case of vanishing pressure (charged dust (CD)) has received considerable attention. The general solution, in which 
the fluid density equals the norm of the invariant charge density, was presented in curvature coordinates by Bonnor 
0. The proof that this relation characterizes regular CD in equilibrium, i.e. in the general static case, was given 
later ^J,|| . In the spherically symmetric case another proof was proposed in Ref . Q . Concrete CD solutions were 
studied in these coordinates MM. The generalization of the incompressible Schwarzschild sphere to the charged 
f***- [ case with constant Tq was also undertaken in a CD environment ||. Charged dust, however, has been investigated 
more frequently in isotropic coordinates, since these encompass the entire static case and allow to search for interior 
solutions to the more general Majumdar-Papapetrou electrovacuum fields |7|,[|. In both coordinate systems there is a 
simple functional relation between 1700 and the electrostatic potential. In isotropic coordinates there is one non-linear 
main equation which has been given several spherical jftJO] and spheroidal [ pl|^2|Jl^ ] solutions. One of them 
coincides with the general static conformally flat CD solution . These CD clouds may be realized in practice by 
a slight ionization of neutral hydrogen, although the necessary equilibrium is rather delicate. They have a number of 
1 , interesting properties: their mass and radius may be arbitrary, very large redshifts are attainable, their exteriors can 
J~ ' ■ be made arbitrarily near to the exterior of an extreme charged black hole. In the spheroidal case the average density 
can be arbitrarily large, while for any given mass the surface area can be arbitrarily small. When the junction radius 



ro shrinks to zero, many of their characteristics remain finite and non-trivial. One can even entertain the idea for a 
point-like classical model of electron were it not for the unrealistic ratio of charge e to mass m jjj] . 

Recently, new static CD solutions were found, in particular with density which is constant or is concentrated on 
thin shells 16 Tl§. In the spherically symmetric case a relation has been established with solutions of the Sine-Gordon 
and the Ai/^equations 18 



The necessary condition for a quadratic Weyl-type relation has been derived also for perfect fluids with non- vanishing 
pressure fjJ]fo|]. However, in this case many other dependencies between the electrostatic and the gravitational 
potential are possible, even when combined with constant Tq 

The original Schwarzschild idea of constant density has been also tested in the charged case for a perfect fluid |2^] 
or for imperfect fluid with two different pressures [[l]. An electromagnetic mass model with vanishing density has 
been proposed in Ref. Unfortunately, the fluid has negative pressure (tension). Although the junction conditions 
do not require the vanishing of the density at the boundary, this is true for gaseous spheres. A model with such type 
of density was proposed both in the uncharged and the charged case [^(|^7| . 

Another idea about the electromagnetic origin of the electron mass maintains that, due to vacuum polarization, 
its interior has the equation of state p + p = 0, where p is the density and p is the pressure. This leads to tension, 
easier junction conditions and realistic e and m [p8H31[, It can be combined with a Weyl-type character of the field 
p2|. The experimental evidence that the electron's diameter is not larger than 10 -16 cm, however, requires that the 
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classical models should contain regions of negative density |33],£54| ■ Probably an interior solution of the Kerr-Newman 
metric is more adequate in this respect. 

The presence of five unknown functions and just three essential field equations allows one to specify the metric and 
solve for the fluid characteristics [p5| . This is impossible in the uncharged case. Another approach is to electrify some 
of the numerous uncharged solutions. This has been done for one of the Kuchowicz solutions |36| in Rcf. [m]]. Two 
other papers [j38|]39f l build upon the Wyman-Adler solution |40|j4l]| . Thus a charged solution is obtained, which has 



approximately linear equation of state when m/ro is small. In Refs. 42-44 generalizations of the Klein- Tolman (KT) 
solution were performed. Recently, static uncharged stars with spatial geometry depending on a parameter 

J47| , |48[ have been generalized to the charged case [Q. General transformation connecting uncharged and charged 
solutions was introduced in Ref. |5(|. 

The purpose of the present paper is to present a new and simple classification scheme for charged static spherically 
symmetric perfect fluid solutions. The calculations in each case are pushed as further as possible and general formulae 
are given in many instances. The known solutions are reviewed and compartmentalized according to this scheme in 
order to illustrate general ideas, without being exhaustive. New solutions are added where appropriate. The intention 
has been to stick to the simplest cases and remain in the realm of elementary if not algebraic functions. The emphasis 
is, however, on the general picture, which appears unexpectedly rich and simpler than in the uncharged case. 

The metric of a static spherically symmetric spacetime in curvature coordinates reads 

ds 2 = e v dt 2 - e x dr 2 -r 2 dil 2 , (1) 

where dil 2 is the metric on the two-sphere and v, A depend on r. The fluid and its gravitation are described by 
five functions depending on the radius: A, v, p,p and the charge function q which measures the charge within radius 
r. There are only three essential field equations, hence, two of the above characteristics must be given. We shall 
classify solutions according to this feature. For example, (y, A) is the case of given metric and the other three fluid 
characteristics are found from the equations. This does not mean that solutions are distributed among groups which 
do not overlap. Thus, (p 7 q) is a completely general case - any solution, after p and q are known, may be put into this 
class. The essence is that p and q 2 are given functions and there is control over them; they can be chosen regular, 
positive and comparatively simple. Then the other three functions are usually more complex and are not always 
physically realistic. 

In Sec. II the Einstein-Maxwell equations are organized into three main and two auxiliary ones. Two of the main 
equations are universal for all cases, the third one varies from case to case. Cases with given v have linear first-order 
differential equations for e~ A . Cases with given A have linear second-order equations for e v l 2 or non-linear first order 
equations for p + p. These results hold also for fluids with a linear equation of state. The difficulties in the cases (p,p) 
and (/?, q) are discussed, which prevent their analytical treatment. The junction conditions are given in general form 
and a reasonable set of physical requirements is included. 

In the following table a summary is given of the discussed cases, the known solutions which are reviewed and the 
new general and particular solutions. They have been checked by computer for realistic properties. 



TABLE I. Summary 



Section 


Cases 


New solutions, Eqs. 


Old solutions, Refs. 


II 


(p,p) , (p,q) 


none 


none 


III 


(A,f),(A,Y) 


36-8 


21,35 


IV 


(A,n=-1) 


41-4, 49, 55 


28,29,31,32 


V 


O.p) 


56-8, 75-8, 81, 84-7 


1,4-6,10-13,15-18, 59 


VI 


0> i) 


92-4 


38,39,63-67 


VII 


{v,n) 


105-6, 108-10, 111-3, 114 


36,53 


VIII 


0,p) 


117-8 


none 


IX 


(A,p) , (A,g) , (p,q) 


136-9, 145-6 


6,22,23,25-27,37,49,64 


X 


(A,n) 


154, 156 


42-44,77 
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The parameter n designates the linear equation of state Eq. (17), which represents a relation between p and p, while 
Y = p + p. Cases which contain q are direct generalizations of uncharged solutions depending usually on a parameter. 
After some reformulation the same property is shared by the cases of charged dust (v,p = 0) and electrified de Sitter 
solutions (A, n = —1). 

Sec. XI contains some discussion and conclusions. 

II. MAIN EQUATIONS AND CLASSIFICATION 

The Einstein-Maxwell equations are written as 

*T ° = K P+ ^ = ^e-* + l(l-e-*), (2) 

"P-$ = T*- X -?(l-e- X )> (3) 



where the prime means a derivative with respect to r and k = 87rG/c 4 . We shall use units where G = c = 1. The 
charge function is obtained by integrating the charge density a. We shall use, however, q as a primary object and 
then 

Ka =% e -V 2 . (5 ) 

When cre A / 2 = const, q <~ r 3 , the case of the so-called constant charge density. Spherical symmetry allows only a 
radial electric field with potential <j> given by 

F 01 =<t>> = -±e^. (6) 



Eq.(2) may be integrated by the introduction of the mass function 

M(r) = ^J\*T°dr = ^J o 

and gives 



This can be rewritten as 



x 2M 

z = e~ A = l . (8 

r 

q 2 2M' 1 , 

*P+^ = -pr = ^(l-z-rz'), (9) 
and constitutes the first of our main equations. The second is obtained as a sum of Eqs. (2) and (3): 

«0, + p) = ^(i/ + A') = V-f (10) 

The third main equation will not be Eq.(4) but another combination of Eqs.(2)-(4) which varies from case to case. 
One can transform Eqs. (2)- (4) into expressions for p, q, p 

_ x (v" AV v n A' 2>v' 1 \ 1 
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r 4 I 2 4 4 2r 2r - 2 / - 2 ' 



\ ( v" AV i/' 2 5A' i/ 1 \ 1 

2«p = e- A -— + — + + +-?• 13 

\ 2 4 4 4r 2r r 2 / r 2 v ; 

These equations may be written as linear first-order equations for z, suitable for the cases (u, q) , (y, p) and (v,p). 
Introducing y = e v l 2 we have 

4o 2 

(r 2 y' + ry) z' = -2 (r 2 y" - ry' - y) z - 2y + -±-y, (14) 

(r 2 y' + 5ry) z' = —2 (r 2 y" — ry' + y) z + 2y — Anpr 2 y, (15) 

(r 2 y' + ry) z' = —2 (r 2 y" + 3ry' + y) z + 2y + Anpr 2 y. (16) 

In the uncharged case the prescription of an equation of state makes the system of field equtions extremely difficult 
to solve. This is true even for the simplest realistic linear equation of state 

p = np~p , (17) 

where n is a parameter taking values in the interval [0, 1] for physically realistic solutions, while po is a positive 
constant, allowing the existence of a boundary of the fluid where p = 0. When po = we obtain the popular 7-law 
(with notation n = 7 — 1). In this case Eqs. (2)-(4) with q = lead to an Abel differential equation of the second kind 
independently from the approach or the coordinate system |5^]. It is soluble in few simple cases. Almost all of the 
few known solutions of the more general Eq. (17) may be obtained by imposing a simple ansatz on A which makes 
the system overdetermined |52|. Therefore, it is surprising that in the more complex charged case fluids, satisfying 
Eq. (17), are subjected to a linear equation, similar to Eqs. (14)-(16). Plugging Eq. (17) into Eq. (10) and replacing 
the resulting expression for p in Eq. (15) yields 

2„ I , 5n + 1 \/_ o/^2// , 3 - n I , \ , o 4K P0 2 



r 2 y' + — ry ) z' = -2 r 2 y" + ——ry' +y)z + 2y f--r 2 y. (18) 

n+1 ) \ n+1 ) n+1 

We call this case (v,n). Eq. (18) was derived in Ref. |33| when po — 0. Eq. (17) is a privileged one, due to its 
linearity. Another realistic equation of state, namely the polytropic one, reads p = np 1+1 / k . It causes the appearance 
of radicals in Eq. (15) and leads to non-integrable equations. 
Eqs. (14)-(16) and (18) may be written in the general form 



gz' = f lZ + f , (19) 



whose quadrature is 



z = e F (C + H), (20) 

^dr, H= [e' F ^dr. (21) 
9 J 9 

Here and in the following C will denote a generic integration constant. They may be written also as linear second-order 
differential equations for y, useful in the cases (A, q) , (A, p) , (A,p) and (A, n): 

2r 2 zy" + (r 2 z' - 2rz) y' + {rz' - 2z + 2 - *^-\ y = 0, (22) 

2r 2 zy" + (rV - 2rz) y' + (5rz' + 2z - 2 + 4/tpr 2 ) y = 0, (23) 

2r 2 zy" + (r 2 z' + 6rz) y' + (rz' + 2z - 2 - Anpr 2 ) y = 0, (24) 
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2r V + ( r 2 z' + 2— -rz] y' + ( ^^-rz' + 2z - 2 + y = 0. (25) 

\ \ n + 1 n + 1 J 

The coefficient before the second derivative is one and the same in all cases. Eq. (22) is the generalization to the 
charged case [^,^9) of the Wyman equation ]io| . Obviously, the case n = — 1 is not covered by Eqs. (18), (25). 

One can find first-order differential equations also for the cases (A,*), based on the well-known Tolman- 
Oppenheimer-Volkoff (TOV) equation |l6|,Q, generalized to the charged case p9| : 

^ 2M + K P r 3 -q 2 /r {q 2 )' 
KP = ' K(/? + P) 2r(r-2M) + ~ ' (26) 

We can trade q in Eq. (26) for p and A by using Eq. (9). The result is a Riccati equation for Y 

Y' = ~\e x Y 2 + -Y -Ak?- + \ (r 2 M')' . (27) 

It marks another way to solve the cases (A, p) , (A, n = — 1). Its solution yields for the pressure up = Y — up. Eq. 
(27) is a non-linear, but first-order companion of Eq. (15). In somewhat different notation it was derived in Ref. fl43| . 
Using Eq. (9) and the definition of Y we find Riccati equations for the cases (A, q), (X,p) and (A,n) with the same 
coefficients before Y' and Y 2 . In the last case it reads 

Y> = - VY 2 + _ 4 \ _ ^Po 2 , 

2 V 2 (n + 1) rj (n + 1) r r 4 V ; V ; 

These equations may be transformed into linear second-order equations by a change of variables from Y to y and we 
obtain exactly Eqs. (22)-(25). In this process we have exchanged Y, which is the sum of the pressure and the density, 
for y which is part of the metric. 

So far we have reformulated the original system of Eqs. (2)-(4) into Eqs. (9), (10) and a third equation, presented 
in many different forms, adapted to the cases of the proposed classification. We have briefly discussed the cases 
(A, v) , (A, *) , {v 1 *). The three remaining cases (p, p) , (p, q) and (p, q) are the most natural ones since one prescribes 
two of the fluid characteristics, hoping that the third one and the metric will be regular and reasonable. The case 
(p, q) is easily reduced to (A, q) because of Eq. (9). In the case (p,p), Y is also known and Eq. (27) becomes an 
intricate, non-linear, second-order equation for M which is not simpler than the TOV equation. It seems that it can 
be dealt with only numerically. 

The third main equation for the case (p, q) is obtained in the following way. Let us replace the density in the TOV 
equation (26) with its expression from Eq. (9). After some tedious manipulations, the following equation for M is 
found 

(M + go) M' = f 2 M + f 3 , (29) 

whose coefficients are functions of p and q. There is a standard procedure for the solution of such equations |^l],[55| . 
It consists of two changes of variables which bring them to the Abel equation of the second kind 

uto c - u - / (C) , (30) 

Its integrable cases are few, depend on the shape of / [Q and are tabulated in Ref. ]5q] . 

As a whole, the most attractive are the mixed cases (v, *) where one fluid characteristic and one metric component 
are specified. 

The five functions which describe the fluid together with its gravitational field should satisfy some physical require- 
ments. Eqs.(7)-(8) show that at the centre M (0) = and e A = 1. The density and pressure should be positive and 
monotonically decreasing towards the boundary. It is obvious that q 2 should be positive, too. The boundary ro of 
the fluid sphere is determined by the relation p (tq) = where a junction to the Reissner-Nordstrom (RN) metric 

x 2m e 2 , , 

e" = e~ A = 1 - — + — (31) 

should be performed. The metric and v' must be continuous there. This leads to the expressions 

^ = l-e^l + ^U^M + ^l I (32) 
r \ 2 J r 2r^ 
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= l-e (l + roi/) = § — 



(33) 



The condition e = q (ro) follows from the vanishing of the pressure and vice versa. 

Solutions, satisfying the above conditions (which is already rather non-trivial) are called in the following physically 
realistic. Some other, more stringent requirements, are discussed in Refs. [^6]-^8|. The most important of them is 

< {j3^j — 1) which means that the speed of sound is positive and causal. One has control over this characteristic 

in CD and models with linear equation of state. The case (v,q), which is a generalization of uncharged solutions, 
shares their behavior for small q. In other cases the expressions for density and pressure may be so complicated that 
the fulfilment of this requirement is hard to estimate. 



III. THE CASES (A, v) AND CONSTANT T ° 



In the case (A, v) Eqs. (11)-(13) should be used to find q,p and p. This is the simplest case but control over pressure 
and density is completely lost and one must proceed by trial and error. Krori and Barua |35fl have given the solution 



A = air 2 



v = a 2 r 



a>3, 



(34) 



where <Zj denote generic constants of a known solution. They are fixed by the junction conditions. The solution is 
non-singular and the positivity conditions are satisfied. 

Eqs.(7)-(8) show that the generalized Schwarzschild condition Tg = const determines A 



1 



or 



(35) 



where a is a positive constant. We have the freedom to choose one more function. In Ref. [pi the condition p = was 
further imposed. We shall review this solution in Sec.V. In Ref. Elj] some relations between cf> and v were utilized. 
They can be arbitrary in the perfect fluid case, so we do not base the classification on such criteria. One should use 
instead Eqs. (22)-(25) or Eq. (27). The simplest expressions are obtained in the case (A, Y) when Eq. (27) is used to 
determine p 



Anp 



-rY' 



r 2 Y 2 arY 



2(l-ar 2 ) 1-ar 2 



12a. 



Then up = Y — up and 



r 4 (3a — up) 



Y' 



rY (Y - 2a) 
2 (1-ar 2 ) 



The function v is determined by a simple integration from Eq. (10) 

v' = —^{Y-2a). 

1 — ar z 



(36) 



(37) 



(38) 



This is an advantage over the second-order equations (22)-(25). Eqs. (35)-(38) solve the problem in a minimal algebraic 
way. Several positivity conditions follow for the master function Y and for p: 3a > up, Y > 2a, Y > up, Y' < 0. 

In order to illustrate how the above scheme works one may take any of the solutions elaborated upon in Ref. [£l| , 
extract Y and check how the above equations and inequalities are satisfied. The simplest case has 



Y = 2a + oi (1 - ar' 



iV2 



(39) 



Of course, in any of the cases (A, *) we can study the subcase given by Eq. (35). This will be done in the following 
sections. 

A similar problem arises when v is given. Which function should be prescribed in addition in order to have the 
simplest algorithm for generating solutions? Again, Y is the best choice. This is seen by taking the difference of Eqs. 
(23), (24) or reformulating Eq. (10) 



yz — 2y'z + ryY = 0. 
This equation is much simpler than any of Eqs. (14)-(16). 



(40) 



G 



IV. CHARGED DE SITTER SOLUTIONS AND THEIR GENERALIZATION 



The uncharged case is a special integrable case with linear equation of state (17) and n = — 1, po = 0, which is 
equivalent to the de Sitter solution. The charged case is also completely solvable. The third main equation is Eq. 
(27) which becomes 

(r 2 M')' 1 

«P=^3^ = ^(2M' + rM»). (41) 
Eq. (10) gives v = — A This is a feature also of the exterior RN solution. Eq. (9) yields 

q 2 = ~(jf)' = l(2M>-rM»). (42) 
Thus, when M is given, all other unknowns follow from simple formulae. Eq. (42) allows also one to take q as a basis: 

2 pr 2 

K p = 2a -^-A i_dr, (43) 

T Jo r 



M = -2 f f 2 f ^tffdf + % 3 . (44) 
Jo Jo r 3 

Here ao is some positive constant and clearly q — r 2+e qo (r) where e > and go (0) = const. Eqs. (43)-(44) 
demonstrate the process of 'electrification' of de Sitter space. The bigger the charge function, the lower the density 
until some point ro is reached where p (ro) = and consequently the fluid sphere acquires a boundary. We already 
know that M > 0. Eqs. (41)-(42) show that p and q 2 are both positive only when M' > and 2M' > r \M"\. The 
equality holds at the boundary, where M" (ro) < 0. The charge and the mass of the solution follow from the junction 
conditions (32), (33) 

e 2 = 2r 2 M'{r ), m = M (r ) + r M' (r ) . (45) 

Obviously, there is an abundance of solutions since M and its first two derivatives have to satisfy few simple 
inequalities. Four solutions are known in the literature. The T® = const condition leads to q — in the interior 
and consequently to the de Sitter solution, which has constant density. One can introduce, however, a surface charge 
a ~ S (r — ro) which gives q — 6 (r — ro) e, p — po and 

z = l- ^p Q r 2 . (46) 

This is the solution of Cohen and Cohen . Another solution |2^] has 

M = ^a V (5a? - 2r 2 ) , (47) 

Here tro is the constant charge density, while a\ is related to ao = t^OqOi. This is one of the electromagnetic models 
of the electron. When ao — > the limit is flat spacetime. We have shown, however, that ao is not obliged to vanish 
when q = so that, in general, the case (A, n — —1) is an electric generalization of de Sitter spacetime. The natural 
generalization of flat spacetime is the CD solution, discussed in the next section. Another solution with M similar 
to the one in Eq. (47) and satisfying the relation np = q 2 /r 4 was given in Ref. It has zero total charge and a 
simple expression for coming from Eq. (3). 

A fourth solution has been found by Gautreau |32j. It has 

2A/ = c 2 r-^isin 2 -, (48) 
r a 2 

and q = —r 2 4>'. Many other solutions are possible and we give one simple realistic example, a variation of the solution 
in Ref. @: 

M = a (r 3 - r 4 ) , 2 = 1- 2ar 2 + 2ar 3 (49) 
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Kp = 2a (3 - 5r) , g 2 = 2or 5 (50) 

The junction at r = 3/5 gives m = ar^, e 2 = 2ar$. 

The case Y — is easy because Eq. (27) collapses into the simple relation (41). A generalization can be made by 
taking Eq.(41) as a basis. Then q is given again by Eq. (42), while Eq. (27) becomes 

Y' = - r -e"Y 2 + ^Y (51) 

This is a Bernoulli equation and, unlike the Riccati equation, it is readily soluble in quadratures. Its general solution 
gives an expression for the pressure 

* P= C + //e3V2 rdr -*ft < 52 > 

where C _1 = Y (0). When C — > oo we return to the previous case Y = 0. Fortunately, Eq. (10) can be integrated 
explicitly too and a closed expression is obtained for v 

e v ' 2 = A- X e- X < 2 (l + -L f e 3X ' 2 rdr\ , (53) 



A = 1 + / ° e 3A / 2 rdr. (54) 
2C Jo 

The second equation follows from the junction conditions. Thus, every function M leads to two solutions: one with 
trivial Y and one with non-trivial Y, satisfying Eq. (51). The trivial solution has the disadvantage of negative 
pressure. Eq.(52) suggests that in the non-trivial case solutions with positive pressure may exist. 
Let us take Eq. (47) with somewhat different constants 

M = br 3 (1 - br 2 ) . (55) 

where b > 0. When x = r 2 < 2/56 the density is positive and decreasing. The pressure is positive at the centre when 
1/46 < C < 5/126. Let us choose 6 = 0.01 and C = 40. Then p has a maximum at xq — 0.3 and a root at xq = 2.5. 
This is a semi- realistic interior solution with e 2 = 56 2 ieq and mass given by Eq. (45). 



V. THE CASE OF GIVEN PRESSURE. CHARGED DUST 



In this section we discuss the cases (v,p) and briefly (A,p). The pressure is considered a known positive function 
which decreases monotonically outwards and vanishes at the boundary of the fluid sphere. The simplest case is p = 0. 
This represents charged dust. The case (v,p) is much easier since the third main equation (16) is linear and first-order. 
Something more, in Eq. (19) /i = —2g'. This allows to obtain a compact expression for z 

1 C 4k fW rv'\ „ , 

z = 2 + -2 + -2 / 1 + -tt) « i r 3 pdr. (56) 

The knowledge of v allows to satisfy two of the junction conditions by choosing a function continuous at rg together 
with its derivative. Eqs.(9)-(10) provide an expression for q in terms of the known v,p and with z (f,p) from Eq. (56) 

a 2 

^ = 1 - z(l +rv') + npr 2 . (57) 

The density follows from Eqs. (10) and (56) 

_ 2z 2 2zv' zv' zv" A K p 

When r — > 0, z — > 1 and the third term on the right produces a pole unless v 1 — ► 2vqt, which means v" — » 2vq. 
Then rv' — > and the poles in the first two terms cancel. The last two terms approach negative constants. In order 
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to compensate them and have a positive density at the centre, vq must be a positive constant and the inequality 
8^0 > & K P (0) should hold. Another consequence is that e v is an increasing function in the vicinity of the centre. 
Since e~ A is a decreasing function, which starts from 1 and equals e v at ro, we have e v (0) < 1. Eq. (56) shows that 
z has a pole at r = unless C = 0. 

Now we can understand the physical meaning of the terms in z. The third term represents the contribution from 
the pressure to the metric. When p = 0, z is still non-trivial and represents the general CD solution. It has a pole at 
the beginning of the coordinates and should not be used there. When C = only the first term remains and this is the 
regular CD solution in curvature coordinates jy. The intricate proofs that this is the most general regular solution 
are replaced here by the obvious condition C = 0. The second term in Eq. (56) may be induced in principle by 
the third if the lower limit of the integral is changed. The first term may be absorbed by the third if one makes the 
shift up — ► Kp + ^2 ■ 

Let us discuss in detail first the case of regular CD. It provides an excellent illustration of the classification scheme. 
In the case (u,p — 0) Eqs. (5)-(6) and (56)-(58) give 

i / rv'\ 2 r 2 v' , . 

V 2 J 2(1 + ^)' y ' 



rv" + 2v> + \rv' 2 ± j 
k P =kg= \ , = -e u/2 + cj) . (60) 

r(l + ^) 3 



For the case (X,p — 0) we should use Eq. (59) to express everything via A 

, 2 



v = — 

r 



(e A / 2 -l), q = r(l-e- x ' 2 ), (61) 



2 

Kp = KO = 



-|(e-V»_ e -A + :A'e-*). (62) 



These formulae involve only the first derivative of A. The case (q,p = 0) is explicitly solvable, unlike the general case 



,A_ r 2 n/f _^ q 2 



e =iTtf' M = q ^ ^ 



v' = r , Kp = na= ^- (r - q) q 1 . (64) 

r (r — q) r a 

The last formulae clearly demonstrate the electrification of flat spacetime, which is the trivial dust solution in the 
uncharged case. If we put p = in Eq. (29) we still obtain a complicated Abel equation. One can check that M, 
given by Eq. (63), satisfies it. The most complicated case is (p,p — 0). Introducing w via e A = w~ 2 , Eq. (61) gives 

/ 2(1 -to) , 
is' = — '-, q = r(l-w), (65) 



rw 



while Eq. (62) becomes an equation for w 



rww' = —w 2 + w — —r' 1 p. (66) 

This is an Abel differential equation like Eq. (29). The procedure for its solution was described briefly after Eq. (29) 
and brings it to its canonical form [5l]j55] |. The set of density profiles, leading to integrable w is very restricted. 

Several explicit dust solutions in curvature coordinates have been given. Efinger || studied the simple case e~ A = 
4/9 in the presence of a cosmological constant. When it is zero we have q = r/3 and singular e v = air, p = a = 4/9nr 2 , 
Fw = \\[&\'r~ 1 l 2 . Florides illustrated his general discussion Q] with a power- law solution q = a2r a3+3 . Finally, in 
Ref. § the case of constant Tq was discussed. It leads to Eq. (35) and its replacement in Eqs. (61)-(62). 

As was mentioned in the introduction, most work on CD has been done in isotropic coordinates 
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ds 2 = U- 2 dt 2 - U 2 (dr 2 + rW) , (67) 
where the simple relation <j> = ztU^ 1 holds and the only equation to be solved is 

U" + -U'= -^U 3 p. (68) 
r 2 

When U is given the density is readily determined from Eq. (68). The function U~ 2 = air 2 + a 2 was used in Refs. 
[[unpin . Another, more general function 

r + kr k +1 [ ' 

was studied for fc = 2 in Ref. [jl2|, for fc = 4 in Ref. jyj and for general fc in Ref. Recently, the function 

U = ai s ^L was studied 

When p is given, it is convenient to transform Eq. (68) into 

U xx = -f^, (70) 

where X = l/r. We may consider p(U) as a known function, chosen to simplify this equation. One possibility, 
leading to Bessel functions, is p — ag/U 2 (l^Jl^]. Other choices lead to two-dimensional integrable models, like the 
sine-Gordon equation |1S| ]. The metric (67) describes also the regular static CD with no symmetry. Then the case of 
constant fluid density po is soluble in cn, one of the Jacobi elliptic functions, but the solution cannot be spherically 
symmetric. In the spherically symmetric case Eq. (70) becomes an Emden-Fowler equation, whose integrable cases 
are tabulated in Ref. [p5| . Eq. (70) is not among them. 

Let us go back to curvature coordinates (1) which are more convenient when one wants to study also the case of 
non- vanishing pressure. From the considerations about the regularity of p, given after Eq. (58), it follows that the 
simplest function v would be v = v$r 2 . This choice coincides with Eq. (34) and leads to exponential functions and 
probably to the error function in Eq. (56). We shall choose a different option 

e»=(a + br 2 ) k , (71) 

where fc is an integer, < a < 1 and b > 0. This ansatz has been thoroughly studied in the uncharged case p6| , p8| . 
In this section we take k = 1. The CD solution build on it has 

+ br 2 \ 2 _ br 3 

+ 2^2/ ' q ~ a + 2br 2 ' { ' 



2b (a + br 2 ) (3a + 2br 2 ) 
k P = — i LA L, (73) 

(a + 2br 2 ) 3 



Density is monotonically decreasing. The junction conditions give 

bv 3 

e=^ r m = r (l-a-2^). (74) 

The condition z = e v at ro supplies the relation Tq = (l — 4a + \/\ + 8a) /8b, which ensures the obligatory e = m. 
The r.h.s. is positive when a < 1, which is the case. 

The main disadvantage of regular CD solutions is that they possess a fixed ratio e/m which is unrealistic, especially 
for classical electron models, and requires the extreme RN solution as an exterior. Eq. (56) tells that when the point 
r = is excluded, general CD solutions are possible. They have the following characteristics 



-2 



-V 



2 \ r 



1 + C-), (75) 



rv \ I „e \ v z 



^=( 1 + - (1 + ^)7-7. (76) 
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rV 2 - 4C (1 + rv') e~»] 1/2 



2r (r 2 + Ce^) 



^1/2 



(78) 



When C = they reduce to Eqs. (59)- (60). What is the physical significance of such solutions? Let us consider a 
core of perfect fluid, occupying a ball with radius ro, with given v and p. Then z is determined from Eq. (56) with 
C 0: 



1 



1 



C(ry 



2 



(79) 



C (r) =4k I [ 1 + — ) e"r 3 pdr. (80) 



o 

At the junction z is equivalent to Eq. (75) where C = Co = C(r ). An observer cannot understand whether the 
interior solution consists of perfect fluid or general CD since the only imprint left by the pressure distribution inside 
is a constant. In principle, the RN metric should be taken as an exterior, but there are no obstacles to take a general 
CD metric and postpone the junction to RN till another point rj > r . Thus we obtain a triple-layered model with a 
perfect fluid core up to ro (zone I) , a halo of general CD from ro to r\ (zone II) and the RN solution for r > r\ (zone 
III). The metric component z is given correspondingly by Eq. (79), Eq. (75) with Co and Eq. (31). For this purpose 
we choose a continuous, together with its derivative, function v in the region < r < r%. At the first junction z is 
also continuous, while the pressure drops to zero. The density is finite. In zone II the pressure remains zero, while the 
density continues to decrease. Finally, at r% the density also drops to zero (perhaps with a jump) and a RN solution 
follows till infinity, making the composite solution asymptotically flat. As seen from Eqs. (76)- (77), the constant Co 
shifts e/m from 1, like a perfect fluid solution, occupying zones I and II would do. 

This picture can be backed by an explicit example with the ansatz (71). Since the details should not depend on the 
form of p (r) but only on the value of C, we shall discuss first zone II and its junction with zone III. In the interval 
[ro, r\] we have 

(a + br 2 ) [C + r 2 (a + br 2 )] 
r 2 (a + br 2 ) 2 

The total mass and charge are given by 

m = n (l -a- 2br\) , e 2 = r\ (l - a - 3brf) . (82) 
The junction condition at r\ gives 

C =m 2 -e 2 . (83) 

This relation clearly shows the effect of pressure in zone I on the mass-charge ratio, not altered by the CD halo in 
zone II. The result in zone I was found in Ref. 

An example of a physically realistic solution is given by a = 0.01, r*i = 4.013v / Co, ro — 3.588%/Co, b — 0.075/VCo- 
This proves the existence of CD solutions with |e| jm ^ 1. They sustain the value of C, obtained in zone I from a 
perfect fluid solution with positive pressure. In Ref. ]59[ ] the same question was studied in isotropic coordinates but 
the solutions have negative density in the vicinity of the origin. 

It seems that there are no solutions of the type (v,p > 0) in the literature. It must be stressed that this case is 
completely general, unlike the cases with constant Tq or with p + p = 0, discussed in the previous sections. Every 
perfect fluid solution may be reformulated as a {v,p) case. Let us proceed with the ansatz (71). We base the discussion 
on C (r) as a fundamental object, hence, it is convenient to select a more involved p leading to a simple C (r) given 
by Eq. (80). Let us take 

KV = (84) 
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and then 



C(r) 



Pa 



(85) 



(a + br 2 ) 
(a + 2br 2 f 



(po 



hir 2 ) 



br 2 



(86) 



(a + 2br 2 Y 



b - bp Q 



-,api 



(87) 



where z follows from Eq. (79) while q is calculated from the general expression (57). Here po and pi are positive 
constants, connected by p = p\r\. The ansatz (71) yields vq = b/a and the condition p(0) > becomes 86/5a > 



up (0) = po/a. This is weaker than a necessary condition for positive q , 6 > po. At the boundary the density is 
always positive. We have also po — 3Co/?"o and using the numerical values from the CD example, one can express 
Po and Pi in terms of Co, in addition to b, tq and r\. It can be shown that for C > 0.06, p and q 2 are positive and 
the solution is physically realistic. Thus we have constructed explicitly a triple solution of the type PF-CD-RN. This 
ends the discussion of case [v,p). 

The case (A,p) relies on Eq. (24). Even the simplest constant z leads to special functions or non-integrable equations 
for y when p is chosen physically realistic. 



VI. THE CASE (v, Q) 

The cases {v 1 q) and (A, q) comprise the most direct generalizations of the numerous uncharged solutions. Setting 
q = one obtains one of them with the chosen ansatz for v or A. Like in the previous section, the case (v, q) is much 
simpler. The third main equation is Eq. (14). Now /i can be represented as a linear combination of p', g and y. This 
is a general result which holds for any of Eqs. (14)-(16), (18). In this particular case we have 

„ ; 8o f v'dr 

h = -2g'+^--Ay, h= ' 



r * J 1 + ^ 



e = — 772 , (89) 



rv' 
2 



H = 2 / e v I 1 + ^- ) e- 2 V- 3 ( % - 1 ) dr, (!>0) 



2 J \r 

and z is given by Eq. (20). We shall use the ansatz (71) where k > 1. In the uncharged case fe = 1 leads to Tolman's 
solution IV pq |. The case k — 2 was discussed first by Wyman |4(| as an extension of Tolman's solution VI and 
later was studied in detail by Adler | fii|. Solution with k = 3 was given by Heintzmann |60) (see also Ref. |5£|). The 
general class was studied by Korkina J61|, but her only explicit solution was the one of Heintzmann. Later Durgapal 
[ p2fl studied in detail the cases k = 1 — 5. All of them satisfy the physical criteria used in this paper, but some have 

1 /2 

irregular behaviour of the speed of sound (dp/ dp) 1 and of p/p. 
Going to the charged case we have 

2 2 / \ m 1 + (fe + 1) txq ( e 2 



e z = q' (so) , - = 1 + — j— f — " 1 , (91) 

r 1 + (2k + 1) txq \Xo 



where r = b/a, x = r 2 . Let us consider first the case k = 1. Then 

1 + ra; q 2 



1 + 2rx (l + C, + 2,Q), Q = £ x3 



dx. (92) 
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The pressure and the density are given by 



t-3C(1 + tx) 2t{\ + Cx) 2{1+tx) f 2q 2 \ AtxQ 

Kp= 1 + 2™ + {l + 2TX f - TT27^ [ 3Q + ^) + (IT^f (93) 



t + C + 3tCx + 2(1 + 3tx) Q q 2 

KP= : ^— " — + +7. 94) 

1 + 2tx x 2 v ; 

When q = these expressions coincide with the pressure and the density from Refs. p6|j6^ ]. We have Kp(0) = 
3 (t — C). Let us choose q 2 = K 2 x 3 . Then the condition p (xq) = yields a negative expression for C which means 
that p (0) is positive. The condition z — e v at xq defines a. Eq. (91) yields e 2 = K 2 x\ and we can express c,a,e,m 
through Xq, t, K. 

Let us consider next the case k = 2. Then Eqs. (89)-(90) give 

z=l + Ce F + e F f 2(a + fa)g2& e F = - (95) 

J (a + 3bx) 1/3 x 3 ' (a + 36a;) 2/3 

When g = this is the result of Adler who used the simplifying assumption e F H = 1 instead of the ansatz (71). 
When q 2 = K 2 x 3 we obtain 

Cx 4a „ , 2 9 9 . 

(a + 3bx) 2/3 56 5 

This is the corrected result of Nduka |3£| who wrote the field equations with an error in the sign of q 2 . In all his 
results K 2 should be replaced by — K 2 and the conclusions correspondingly altered. 

The ansatz (71) with r 2 replaced by r or r 3 was also studied p3| . In another paper Singh and Yadav J39| simplified 
Eqs. (14), (19) by demanding fi = —fa. Letting q = Kr they obtain the Euler equation. It has three types of 
solutions, all of which are singular at r = 0. If K = the non-singular solution of Adler is again recovered. 

Other solutions with singular y were discussed in Refs. |]6^-|66|. In Ref. a charged version (with constant charge 
density) of Tolman's solution V was obtained. Further examples of (i>, q) solutions are given in Refs. |36|,^37j but their 
properties are not examined. 



VII. THE CASE (v,N) 

This is the case with given v and the perfect fluid satisfying the linear equation of state Eq. (17). Realistic solutions 
have < n < 1 which means constant speed of sound and causal behavior. Surprisingly, the charged case is much 
easier than the uncharged one. The third main equation is Eq. (18) with -1. The functions F and H read 



r -(a+l) J3ht 



(97) 



// =•_>/,' I A + r a e~^ (l - ) ,/,, ! !).x i 



where 



f v'dr . 5n + 1 . „. 

hi = / r, A = , 99) 

J A+?f' n + 1 ' V ' 

l-3n „ 4n(9n+l) . . 

a= , [3 =- \ , ' - . (100 

5n + l (n+1) (5ra+l) 

We shall study the range < n < 00 and some negative n. The case n = — 1/5 demands special treatment commented 
upon at the end of the section. The case n = — 1 was discussed in Sec. IV. Each of the coefficients a, j3 and A has the 
same limits when n — > ±00. When < n < 00 their ranges are —3/5 < a < 1, 1<A<5, 0</3< 36/5. Expressions 
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(97), (98) are generic for the cases (v,q) , {v,n) , {v,p). The case (u, q), discussed in the previous section, is obtained 
by putting a = — 3, f3 = 2, A = 1, changing the sign of H and an obvious replacement for q instead of po. The case 
(v, p) will be studied in the next section. 

In the present case Eq. (100) shows that a + 1 is always positive for n 6 [0, oo]. Therefore, z will have a pole unless 
C = 0. Hence z = e F H. After z is found, Eqs. (10), (17) allow to extract the density from the metric 



1 



Kp = 



npo + - \v 

r V z 



(101) 



It is more convenient to express p and q from Eqs. (17), (57) respectively and lay all the difficulty in the calculations 
on the density. Plugging Eqs. (97), (98) into Eq. (101) we obtain 



(n + 1) Kp = npo + 



4npo 



z {rv')' 



z 
A 



1 



(A + rf)' 



A + ^ 




(102) 



which is an analog of Eq. (58) and leads to similar conclusions for v. The poles in the two terms in the square brackets 
cancel when r — > because z — > 1 and v' — > 2vqt. This condition makes p a well-defined function. We shall try for v 
an ansatz which generalizes the one in Eq. (71) 



e v = (a + br s f , 



(103) 



where s > 2 is not necessarily an integer. When s > 2 it can be shown that at the centre the relation p (0) + 3p (0) = 
holds. This is exactly the equation of state for the Einstein static universe (ESU) [^,^8). The pressure is negative at 
the centre, which is unrealistic. Therefore we set s = 2. The integral for z is simplified considerably when 



k + 5 + 4y/k (2k + 1) 



31fc - 25 



(104) 



The first few values are ni ~ 2.15, ni w 0.53, f=a 0.39, = 1/3. The first is beyond the physical range, the fourth 
is an exact number. One cannot probe the region n < k, 0.21 with this method. A series of models is obtained, 
parameterized by k. They have 



/(i + tf-^-ifi-^U 



t^r [l + tf- 2 [A 



[A + k) t] 



(105) 



where t = rr 2 , p = Kp /r. Eq. (102) becomes 



(n + 1) — = p 

T 



ip (1 + t) 



(n+l)[A+(A + k)t] 



2kz 
T+t 







l+t 



l+t 



A + (A + k) t 



L A + (A + k)t 

Akz 

(1 + t) [A + (A + k) t] ' 



(106) 



Now additional terms contribute to p (0) and p (0) and they can be made positive by a proper choice of r and pq. The 
pressure is given by Eq. (17) and the charge function by Eq. (57) which becomes 



rq 



l + (2fc + l)i Kp 
z H 1 

l+t T 



(107) 



An interesting case is k — 4 with — 1/3 exactly. When po = and q = this is the equation of state of pure 
incoherent radiation and the case is not integrable [5l| . This specific charged version, which is not the only possible 
one, however, is integrable and assuming p = 2/3 we have 



1 + ff 



2 .3 



1+4 



(l + t)-(l + 3t) 



(108) 
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up _ 5 + 9* 3(3 + 7*)z |+3*+f* 2 + ±t 3 2 

~ ~ 12 (1 + 3t) + 2 (1 + t) (1 + 3t) ~ 4 (1 + t) 2 (1 + 3t) ~ 3 ' ( ] 

^..(l^. + S,). ( no, 

The pressure is positive, monotonically decreases and vanishes at to = 0.54. The charge density has a negative part. 
The junction conditions yield a = 0.3, b — 0.16/rQ, m = 0.5ro, e = 0.45ro. We have e/m < 1. Probably, all models 
in the series have a range of fi where they are physically realistic. 

The method used to simplify Eq. (98) resembles the Korkina-Durgapal method, which was extended to the charged 
case {v, q) in the previous section. There are other ways to simplify the integral for z even when k = 1. One of them 
is to take again n = 1/3. Then a = 0, A = 2, (3 = 3/2, 7 = 1/2. We have 



l+t 



t 1 / 2 (2 + 3t) 



3/2 



- j (m - 4 + 3/xt) V(2 + 3i)t + -j^V31n (3* + 1 + v^t (2 + 3i)) 



(111) 



where the integration constant was chosen so that z (0) = 1. The pressure and the charge are given by Eq. (17), (102) 
and (107) 

up 3(1 + 2*) 9z 1 ( z l+*\ 



r 4(2 + 3t+ 4(2 + 3*) 2* \2 2 + 3* 

The pressure is well-behaved for any /j, and vanishes at to. However, q 2 also has a zero, but at *i, and becomes 
negative in a region where t > t\. Unfortunately, for most values of /i the inequality t\ < to holds, meaning that 
q 2 < in some part of the interior. This is grossly unrealistic. We have been unable to find by computer simulations 
any realistic solution. 

Another case with k = 1 which leads to simple functions can be guessed from Eq. (100). When n = —1/9, (3 = 0. 
We have also a = 3, A = 1/2, 7 = 1. Although the pressure is negative, it is worth being discussed. We have 



l + t 
(1 + 3*) 2 



3 \ 27 2 



(114) 



When < /i < 4/27 the pressure monotonically increases, crosses the * axis and forms a boundary of fluid with 
tension inside. Let us take /1 = 0. The charge function q 2 is positive and increasing. The junction conditions give 
a = 0.55, b = 0.11/r-Q, m = 0.23r , e = 0.35r where t = 0.2. This solution resembles the classical models of the 
electron in Sec. IV. 

Finally, some words should be said about the case n = —1/5, which is really very special. Let us choose the simple 
ansatz (71) for v. The case is formally soluble but z has poles at the centre and is ill-defined. When po = and 
(7 = 0, this case is soluble |3l[] . It is connected by a Buchdahl transformation to the de Sitter solution with n = — 1 
p8| , p9| and is best expressed in isotropic coordinates. This transformation, however, does not work for electrostatic 
fields and probably is not applicable to the present case of a perfect fluid plus an electrostatic field. 

There are two papers which come near to the results of the present section. A solution with singular y, p and p was 
found in Ref. J53| when po = 0. Another one is contained in Ref. J38| which belongs to the (u, q) case and is mentioned 
in Sec. VI. When m/ro is small, the solution approximately satisfies Eq. (17) with n = 3/10, a value which is rather 
close to 1/3. Its connection to our solutions is not clear. 

A final remark is in order. It may be tempting to take the limit e — > in one of the solutions obtained here and 
try to find analytic expressions for the non-integrable at present cases of 7-law uncharged solutions |pl[ . However, for 
the ansatz (71) used in this section, this leads to ro — > and to flat spacetime. In the case of general i/, it seems that 
the above limit will produce charged solutions with zero total charge. This is suggested by the behavior of tr in some 
of the models, namely, its change of sign in the interior. The proper non-charged limit seems to be q (r) = 0. Then, 
however, we obtain the field equations of the uncharged fluid. Hence, nothing is gained in comparison with the three 
approaches described in Ref. They lead to Abel differential equations of the second kind with few integrable 



15 



VIII. THE CASE (y, p) 



This case arises as a special limit of (y, n) when n — > oo, po ~~ * °°- One sees from Eq. (17) that in this limit 
p — ► po = Po/ n which can be made a finite constant. More generally, the third main equation (15) is obtained from 
Eq. (18) plus the change ^ p{r). Eqs. (99)-(100) give a = -3/5, A = 5, = 36/5. The constant C in Eq. (20) 
is zero again in order to prevent a pole at the centre. Eqs. (97)-(99) 

2ef * f „ /„ n/ \ r „ 3/5 



r 2/5 e . (5+^) 2 



^ r -3/5 e -f »3 _ 2K(0r 2) ^ (n5) 



(116) 



2 



The parameters n and po have been exchanged effectively for the density function and in this respect Eq. (115) 
resembles Eqs. (89)-(90) from the {y, q) case. When p is constant we obtain a generalization of the Schwarzschild 
interior solution (70|. It is much more complicated than the other generalization with constant Tg . As seen from Eq. 
(115), the case p = is non-trivial and leads to purely electromagnetic mass models. 

Let us try again the ansatz (103) for v. The method which lead to an infinite series of realistic solutions in the 
{y, n) case is rather useless here even when s — 2. Therefore, we put k = 1. Then 7 = —1/5 and 

1 K (117) 



£!/5 (5 + 6t) 4/5 7 £4/5(5 + 6t) 1/5 ' 
When p = po the integrals bring forth the hypergeometric function 

1 + f [ (1 1 6 6A kpq (\ 6 11 U 

Z ~ (1 + | t ) 4/5 [ U'5'5 ; ~yJ U'5'T ;_ y 

The charge is given by Eq. (113), while the pressure reads 



(118) 



up 2z up , 

— = — -2z t -^. 119 
r 1 + 1 t 

Searching for an elementary solution, one may consider some sophisticated density profile, leading to a simple 
integral in Eq. (117). Unfortunately we have found only solutions with q 2 < 0. 



IX. THE CASES (p,Q), (A,p) AND (A, Q) 

Having in mind the results in the previous section, it is no wonder that the case (p, q), where there is control over 
the evasive charge function, is much more popular in the literature. Eq. (9) easily relates this case to (A, p) and 
(A, q), and Eq. (22) or (23) may be used as a third main equation. Like {v,q), these cases are direct generalizations 
of uncharged cases. Eqs. (7), (8) define z in terms of p and q 

Z = l-~ f pr 2 dr - - f irdr. (120) 
r Jo r J r z 

Eq. (9) shows that p is regular when z — 1 ~ r 2 for small r. We accept the ansatz 

z = 1 - ar 2 - br 4 . (121) 

It follows from Eq. (120) when q — Kr l with I = 2, 3 and p = po — piT 2 , where po, pi are non-negative constants. The 
ansatz for p goes back to Tolman's solution VII. When I = 2 or / = 3 we have respectively 

z = \-\{npo + K 2 )x + ^x 2 , (122) 
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2 = 1 — + ~ - K 2 ) x 2 . (123) 

The total charge and mass are given by 

e = Kx 3 /2 , (124) 

m = i [a+ (b + K 2 )x ] xl' 2 . (125) 

When pi =^ this is a model of a gaseous sphere, both density and pressure vanish at the boundary |p6| , p7f. We have 
Po = PiXo- When p\ = 0, po ^ this is a generalization of the incompressible fluid sphere for I — 3 |22[ | or I = 2 
pUfl . Finally, when p\ = po = we come to a model with electromagnetic mass ||. The solutions given by Eqs. 
(46)-(47) also belong to this ansatz. The case 6 ^ appears to be simpler, due to a peculiarity in Eq. (22) and leads 
to elementary functions. Eq. (35) demonstrates that the other case, 6 = 0, belongs to the family of solutions with 
constant Tg . Eq. (22) becomes in terms of x 



{I -ax- bx 2 ) y xx - (bx + | ) y x - Q + y 



= 0. (126) 



The coefficient before y is constant when I = 3 (or K — 0). Let us study this case first. Eq. (123) gives 

' b - + ^-\ = J- 
4 2a: 3 / 20 



d^-[ b - + 1 Q=±-( K p 1 -llKi). (127) 



A change of variables brings Eq. (126) to 

y« + = 0, (128) 

e =r - dx . d29) 

Jo v 1 — as — 6a; 2 

This integral has several different expressions, but many properties of the solutions depend only on the fact that 
£ (0) = 0, £ and £ x are positive. Eq. (128) is easily solved and the last function to be determined, the pressure, is 
given by Eq. (10). At the junction the pressure vanishes when 

2y € (a; ) = z x (x ) (130) 

is satisfied. The charge is given by Eq. (124) while the mass becomes 

— + (6^ - Kpx) —J x ' . (131) 



1 

m = — 
2 



When pi = the mass is positive. When pi ^ this assertion still holds because then p±xo = po. The positivity of 
the mass does not impose any conditions on the constants po, £0 and K which determine the solution. The pressure 
reads 

up = 4\/i^- - -np + - (3kpi + 2K 2 ) x. (132) 
y 3 5 v ' 

To be more concrete, let us introduce the three classes of solutions of Eq. (128) for d — 0, d < and d > respectively: 

y = C x + C 2 £, (133) 

y = C ie -^ + C 2 e^, (134) 

y = d sin + C 2 cos (135) 
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The conditions y > and p(0) > lead to C\ > 0, C 2 > when d > and to C 2 > 0, C 2 > Ci when d < 0. 
Therefore, e" increases with r until it meets z which decreases. We have y (0) < 1. 

1) Case d = Q. Eq. (127) means npi = IIK 2 , while Eq. (123) gives 6 = -2K 2 . We also have x = n^and the 
solution is determined by two constants, po and K. Eq.(130) results in 

C 2 = tt^kpo, d = -^Kpofo + v^O) ( 136 ) 



^0 = 1 - 77" 



11 Itfl 



(137) 



and y is given by Eq. (133). C\ is real when v < y/3/5. The charge is given by Eq. (124) and the mass is 

4 3/ 

m = -t^kpqXq 



j^Kpox^ 2 , leading to the ratio \e\ jm = 3/4w > \/l5/4. Due to the negative 6 Eq.(129) reads 



1 



Jn 



a + 26a; 
2V^& 



1 



-6 



In 



1 



2y^b 



In Ref. p7| another expression was used, true only in a part of the allowed interval for v. The pressure reads 



(138) 



(139) 



The condition p(0) > is fulfilled when Ci < 2/11. This holds when v varies near \/3/5 = 0.77. Thus, there are 
solutions with positive pressure. 

2) Case d < 0. This means npi < 11K 2 and the limiting case p\ = is possible. We shall discuss first the general 
case. The boundary conditions fix the constants Ci,2- The pressure is given by Eq. (132) 



up 



4V^dz 



1 



1 



; K Po 



(3/tpi + 2K 2 



C 2e v=^ + Cie-v^? 3 5 
A necessary condition for positive pressure at the center is C% > C\, satisfied when 

np 2 



z x (xo) 



[K 2 - Kpx) X 



(140) 



(141) 



is positive, which leads to np\ > 6K 2 . Simplifications are possible when d — b or C\ — but they don't give solutions 
with positive pressure. 

Let us discuss the subcase p\ = 0, i.e. p — pa p5[ |. Now Eq. (141) shows that z x (xq) < 0, therefore, C 2 < Ci and 
the pressure is negative in this subcase, which was not realized in the above reference. Eq. (124) holds as it is, while 
Eq. (125) reads 



1 



3/2 



We have b = K 2 /5 > and consequently a 2 + 46 > 0. Then 

1 . a + 2bx 1 

t = — = arcsm — , ■= ; 

Vb Va 2 + Ab Vb 



Va 2 + 46 



(142) 



(143) 



In Ref. [^7| an incorrect expression was used. Formula (143) allows to make connection with the results of Ref. ||, 
where in addition po = 0, so that the mass arises entirely from the electrostatic field energy. Eq. (141) still shows 
that the pressure is negative. Eqs. (124) and (142) give the ratio |e| jm = 5/3 \K\ Xq. Now a = and z is positive at 
the boundary when |e| jm > v5/3. The second term in Eq. (143) vanishes and combining it with Eq. (134) yields 



y = C'i exp I — \ — arcsin Vbx ] + C 2 exp [ \ — arcsin Vb: 



(144) 



The passage to arccos scales and interchanges their places, giving Eq. (4.6) from Ref. 
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3) Case d > . Now we have Kp± > UK 2 which permits the limit K = 0. The two cases do not differ essentially. 
The value of z x (x ) is positive, while b < 0. Hence, £ is given by Eq. (138), like in Ref. [p6| . The positivity condition 
is even more complicated than in the previous case and simplifications do not seem possible. This is true even for the 
uncharged case, discussed in the mentioned reference. We shall remark only that the simple case np\ = upo = xq = 1, 
K = has positive pressure profile. For small enough K, the charged case must behave the same way. 

Let us study next the case I = 2 when z is given by Eq. (122). Eq. (126) is not soluble unless b = —Kpi/5 = 
meaning that the density is constant. Then it becomes a particular case of the hypergeometric equation 



x (x - 1) Vxx + iivi + m + 1) x - m] v x + mmv = o, 

where x — ax i Vi + V2 + 1 = 1/2, % = 0, 771% = K 2 /2a. These relations give 



11/ 8K 2 
7?1 = -4 ± IV 1 - — 



(145) 



(146) 



and 771,2 € [—1/2,0). The reality of r]i is ensured by Kpo > 23K 2 . Thus, in the charged case it is not possible to 
put p = pa = 0. This case was studied by Wilson |23|| who did not recognize the appearance of the hypergeometric 
function. He assumed that the pressure was positive and developed a series expansion for y. A linear combination 
between the two fundamental solutions of Eq. (145) is a candidate for a regular metric. In this paper the emphasis is 
laid on solutions in elementary functions, so we shall not pursue this issue further, except in the case where K = 0. 
Then Eq. (145) is easily solved 



2(1 



ar 



,1/2 



«1 



This is the expected Schwarzschild interior solution |7(J| , contrary to the claims in Ref. |2J 
A model for a superdense star with the ansatze 



1 - a 2 r 2 
1 + aia 2 r 2 ' 



Ka 2 r 3 



1 + aia 2 r 2 ' 



(147) 



(148) 



has been discussed in the uncharged case for a\ = 2 |7l] |f2|j , a\ = 7 |7J| and a set of discrete values for a\ ]47j. Recently, 
it was studied for arbitrary a\ both in the uncharged [)48| and the charged case |49fl . Introducing the variable 



/ \ 1/2 

/ Ol \ /, 2\l/2 

(l - V 2 ) Vnvv ~ VUnn + d iVv = °> 



one obtains 



(149) 



(150) 



2K 2 

d x = ai +2 , (151) 

ai 

where d\ > is required. 

Formally, Eq. (150) is a subcase of Eq. (126) with a = 0, b = 1, x — > r, y — > y v , so we can use the machinery 
developed there. We have £ = arcsinry. We may change the variable to 8 = arccosry because Eq. (128) is invariant 
under this change. Hence, y v is given by Eq. (135) and y may be found by integration. Use of trigonometric equalities 
helps to find an expression with two terms 

/ cos [(Va\ + l) S + a 3 ] cos [(y/fa - l) 6 + q 3 ] \ 
V ~ aA \ ^a\ + l V3T-1 r 



This formula was found in Refs. [[18|j49j by resorting first to Gegenbauer functions and then to Chebyshev polynomials. 
The peculiarities of Eq. (126) permit a straightforward derivation in elementary functions. 

Let us discuss next a solution of Nduka |37) , belonging to the class (A, q). He chose a constant z ^ 1 and generalized 
the uncharged solution of Ref. (3(| by taking 1 = 1, i.e. q = Kr. Then Eq. (22) turns into the Euler equation and y 
and p are singular at r = 0. 

The ansatz (121) was examined also in Ref. |Q where besides q a generating function G was given, satisfying the 
relation G 2 zq = const. Therefore this is again the case (A,q). Several solutions were found, including the one of 
Nduka |3q] and a singular in y solution, which is a generalization of the solutions in Refs. 
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X. THE CASE (A, N) 



In this case the fluid satisfies the linear equation of state (17) and some ansatz for A is given. The other metric 
component is found from Eq. (25) - a linear second-order equation for y. 

Let us study first the case z = c < 1 which leads to singular pressure and density, but provides the opportunity to 

(25) becomes 



generalize the well-known KT solution W5 46 5lJj7q 



Eq. 



r 2 y" + k ry' + (-fci + k 2 r 2 ) y = Q 



(153) 



3 — n , 1 — c 2«p 

— TP fcl = ' fc2 = 7 — 7TT- 

n + 1 c [n + 1) c 



Its solution is given by Bessel functions 



y = r 2 



CiJ fe3 (Vfer) +C 2 y fc3 (Vfer) 



(154) 



(1-fcof + 4fci 



1/2 



When po = this becomes 



V 



dr 1 



-k +2k 3 



(155) 



It can be checked that 1 — ko + 2k$ > 0. We have included only the non-divergent at r = solution. The density and 
the charge function are singular at the origin. This solution was derived in Ref. jjl2) and rediscovered in Ref. [j?7) . It 
was generalized by Tikekar [[l3| whose solution belongs to the (Y, n) case and is singular at the origin too. Another 
generalization was offered in Ref. |44| but it alters many of its properties, including the equation of state. 

Let us discuss finally the ansatz (35), imposed by the condition of constant Xq, like in Sec. III. Introducing x = ar 2 
turns Eq. (25) into 



x(x — l) y xx + 



2(n + l)' 



y x + Dy = 0, 



(156) 



D=- r L-(to+l-!^ ). (157) 
2 (n + 1) V a J y ' 

This is once again the hypergeometric equation and its solutions depend on the sign of D. Particular cases are the de 
Sitter solution and ESU. There are many elementary solutions but they have either singularities or it is impossible to 
ensure that both p and q 2 are positive. 



XI. DISCUSSION AND CONCLUSIONS 



Charged static perfect fluids have attracted considerably less attention than the uncharged ones, the number of 
papers being roughly an order of magnitude smaller. The introduction hints how diverse were the approaches to 
their study. Charged dust occupies the first place in popularity. The rest is a mixture of generalizations: of the 
Schwarzschild idea about incompressibility, with the limiting case of vanishing density and electromagnetic mass 
models, generalization of the idea about vacuum polarization which brings forth the de Sitter solution, generalization 
simply of well-known uncharged solutions like those of Adler, Kuchowicz, Klein, Mehra, Vaidya-Tikekar and others, 
generalization of Weyl type connections, well studied in the electrovac case. 

In this paper we have tried to show that the charged case has a life of its own when subjected to a natural 
classification scheme. Surprisingly, in many respects it looks simpler than the uncharged case. The presence of the 
charge function serves as a safety valve, which absorbs much of the fine tuning, necessary in the uncharged case. 
The general formulae derived here show that the abundance of solutions is probably bigger than in the traditional 
case. The proposed scheme becomes rather trivial there, representing ansatze mainly for A, v and sometimes p. In 
the charged case, however, it allows to sort out the different ideas mentioned above. Thus constant Tq leads to the 
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cases (A = 1 — ar 2 , *) , while constant density - to the much more difficult cases (v, p) and (A, p), similar to [v, n) and 
(A, n). Electromagnetic mass models are subcases, often spoiled by negative pressure. The point-like idea seems to 
be viable only for CD, leading to flat spacetime when the pressure does not vanish. The generalization of uncharged 
solutions firmly occupies the cases (X,q) , (p,q) and {v,q). Models with n= — 1 arc intimately related to the soluble 
case n / -1 and are much richer than their traditional protagonist - the de Sitter solution. The other models with 
negative pressure, which seem worth being studied, are the generalizations of ESU with n = —1/3. Finally, the Weyl 
type relations, so successful in electrovac and CD environments, seem rather out of place in a 'pressurized' perfect 
fluid. 

Another advantage of this classification is that it delineates the degree of difficulty and the most convenient points 
for attack of the problem. The easiest cases seem to be (A, Y) and (u, Y). If we insist on authentic fluid characteristics 
and not on their combinations, then {v,p) and (y, q) are the best candidates. On the other hand, the most difficult 
are (p,p) and (p,q), which seem to be accessible only numerically. The most unpredictable case is probably (A, v) , 
since there is no control on exactly those characteristics which must satisfy the majority of regularity and positivity 
conditions. 

A different division is between general and special cases. As was noticed in the introduction, a known solution 
belongs to any of the general cases; a solution of (v,q) may be written as (A,p) or even as (p,p) solution. Only the 
simplicity of two of the five functions v, A, p 1 p 1 q betrays where was the starting point. The special cases that we 
discussed, apart from the ansatz (35), are essentially three: p + p = and (A, n), {y,n). All of them have a linear 
equation of state. It should be mentioned that all integrable uncharged solutions with 7-law found in Ref. |[l[] have 
their charged generalizations with p ^ and, hence, possess a boundary. The case n = 0, po = represents CD and 
arises as an electrification of the trivial uncharged dust solution (flat spacetime) which is seen from Eqs. (63)-(64). 
Its boundary may be put anywhere, since p = everywhere. The case n = —1 generalizes the de Sitter solution 
into a bunch of new solutions, characterized by the mass function, subjected to several mild restrictions. It was 
completely solved in Sec. IV. The de Sitter solution appears also as a special case in Sec. X . The option n = —1/3 
in the uncharged case is a mark of ESU. It appears in Sec. VII when b = in Eq. (103). ESU appears also as the 
special class D = in Sec. X. The case n = —1/5 is special in the charged theory too, but its metric is singular as 
was mentioned in Sec. VII. The KT solution is generalized in Eq. (155). In the uncharged case these exhaust the 
integrable cases of a complicated Abel differential equation of the second kind. In the charged case it is replaced by 
a linear first-order Eq. (18), which has many other solutions, discussed in Sec. VII, or by a linear second-order Eq. 
(25), discussed in Sec. X. 

From a mathematical viewpoint, the charged case delivers a surprising variety of equations like those of Euler, 
Bernoulli, Riccati, Emden-Fowler, Abel and the hypergeometric equation. The cases containing A quickly lead to spe- 
cial functions whose spectrum is also rich: Ei, cn, Bessel, hypergeometric, incomplete /3-functions, Jacobi, Gegenbauer 
and Chebyshev polynomials, etc. Therefore, a more thorough study would require computer simulations of special 
functions, which are easier than purely numeric simulations. 

From a physical viewpoint, the most interesting results are, in our opinion, the following: 

1) The case of constant Tq is solved by a simple algorithm involving algebraic operations, one differentiation and 
one simple integration upon a generating function Y , which satisfies few simple inequalities, Eqs. (36)- (38). 

2) The charged de Sitter case is completely soluble in terms of the mass or the charge function, Eqs. (41)-(44). 
There is a generalization with positive pressure given by Eqs. (52)-(54). 

3) The general solution for a given pressure consists of three contributions to z; from regular CD, from general CD 
and from the pressure, see Eq. (56). There is a general CD solution which creates a halo around charged fluid balls 
and may postpone their junction to a RN solution. 

4) The solutions of Korkina-Durgapal possess realistic charge generalizations given by Eqs. (88)- (90). 

5) Fluids with linear equation of state are integrable for any n as seen from Eqs. (18), (25), (28), (97)-(100). 
Elementary solutions may be found by several simplification techniques. Physically realistic is an infinite series of 
models given by Eqs. (104), (105)-(106). The status of the model given by Eq. (Ill) is still unknown. 

6) The case (v, p) is closely connected to (v,n) and is much more difficult than its companion with constant Tq, 
even when the density is constant or zero. 

7) The (p, q) case discussed in Ref. [^7| is incomplete and has errors. It is treated in detail in Sec. IX and its 
connection with Ref. || is elucidated. It has several realistic subcases, but electromagnetic mass models all seem to 
have negative pressure. A formal, but intriguing parallel is drawn to star models with simple spatial geometry 49 1. 

8) Solutions of the case (A, n) with constant Tq are expressed in hypergeometric functions, see Eqs. (156), (157). 
All degenerate elementary solutions have either negative pressure or singularities. 

A more detailed version of the present paper can be found in Ref. |78] . 
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